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Abstract 

We study a discrete-time random walk on the non-negative integers, 
such that when is reached a jump occurs to an arbitrary location k > 
with probability pk- We obtain an asymptotic formula for the expected 
position at time n, in dependence on the probabilities pt and on the 
starting position. Our proof of this result displays the relevance of the 
spectral analysis of the transition operator associated to the stochastic 
process, both of its eigenvalues and of its resonances. 



1 Introduction 

We consider a discrete-time random walk performed on the non-negative integers 
{0, 1,2, • • •} such that, 

• When at position j > 1, one moves left or right with equal probability, 

• When at position j = 0, one jumps to an arbitrary integer k, with proba- 
bility pk , where 

oo 

Pk>o, J2p* = 1 - ( L1 ) 

fc=0 

In other words, denoting by X n the position at time n: 

P(X n+1 =j-l\X n = j) = P(X n+1 = j + 1 I X n = j) = 1, j > 1, 
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P(X n+1 = k I X n = 0) = Pk , fc>0. 



This process is a special case of a more general class of processes considered by 
Lalley [5]. Specialized to the case considered here, the main result of [5] is that, 
under the condition po < i, for each j, k there exists a constant Cjk so that, as 

n — > oo, 

P(X n = k \ X = j) C jk ^=, 
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where ~ means that the quotient of the two sides approaches 1. Drmota ([3], 
Theorem 1) proved a more precise result in the special case in which pk = for 
k > 2. 

Here we are interested in obtaining a more precise characterization of the de- 
pendence of the asymptotics on the parameters defining the process, namely 
the jump probabilities pk and the starting point j. However, we shall not ana- 
lyze the probabilities P(X n = k | Xq = j) themselves, but rather the expected 
position at time n, assuming the starting position Xq = j, 

oo 

E(X n \X =j) = J2 kp (X n = k\X = j). 

fc=0 

We prove the following asymptotic formula, as n — > oo 

e{Xn\x = i) = ^V-n + *Q^ (1.2) 
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1 E(Y 3 ) <E(Y 2 )-1) 2 E(Y 2 )-1 

j t7R7\ l-J 



-3 3E(Y) 2(E(Y)) 2 J E(Y) 
where Y is a random variable expressing the size of a jump, that is 
E(Y = k)= Pk , k = 0,1,2,- •• 

so that 

oo 

E(Y m ) = Y / Pkk m - 

k=0 

Several features of the above asymptotic formula are noteworthy: 
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• The dominant term, of order i/n, reflects the universality in the behavior 
of the process, as it does not depend on the jump probabilities pk, nor on 
the starting position j. 

• The second term, which is a constant shift independent of n, reflects the 
dependence of the process on the jump probabilities pk- This term can be 
either positive or negative depending on whether E(Y 2 ) > 1 or E(Y 2 ) < 1. 
Note that it too is independent of the starting position j. 

• The third term, of order -y^, already depends on both the jumps at the 
origin and on the starting position j. One can thus say that the memory 
of the initial position decays like . 

The performance of formula (|1.2j) is illustrated in Table [TJ where the results 
given by the asymptotic formula (and its truncation to one or two terms) are 
compared with the precise expected values, computed using the generating func- 
tion derived in section [21 

Formula (ll.2j) is proved under a few restrictions on the jump probabilities {pk}, 
some of which are essential, in the sense that the formula is not valid without 
them, and some technical, in that we believe they can be removed by elaborating 
the proof, but shall not do so here. We now discuss these restrictions. 
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n 


I 


I + 11 


I + 11 + III 


E(X n \X =j) 


10 


2.52314 


3.63425 


5.40355 


5.26359 


20 


3.56826 


4.67937 


5.93046 


5.88291 


50 


5.64191 


6.75302 


L54430 


7.53373 


100 


7.97885 


9.08996 


9.64944 


9.64614 


200 


11.28379 


12.39490 


12.79053 


12.78946 


400 


15.95769 


17.06880 


17.34856 


17.34820 



Table 1: Comparison of results from the asymptotic formula of Theorem [T] 
with the exact value of the expected position after n steps, computed using the 
generating function. I, II, III denote the three terms in the asymptotic formula. 
In this example Po — jq,Pi = p-2 = jq,P3 = ^, and j = 5. 



First of all, for (|1.2I) to make sense, we need that E(Y) > 0. Clearly the only 
case in which E{Y) — is when po — 1 and p^ = for k > 1, so that when 
the origin is reached, one remains there at all times. In this exceptional case we 
will show (see section [5]) that the exact expected value is given by the following 
somewhat surprising result: 



Proposition 1 If pa = 1 then 

E(X n \ X =j)= j, Vj,n>0. 



Therefore, to exclude this special case, we make the assumption 
(Ai) po < 1. 

Another family of exceptional cases is when p/. — for all even fc, so that the 
jump at is always to an odd value of k. In this case we will show that the 
formula (jl.2l) is replaced by 



^\X Q =j) = ^ + ^^ (1.3) 
1 F l E(Y*) , {E(Y 2 )-lf , E(Y*)-1 | j2 | 1 1)n+ , +1 i 1 



O 



2tt 
1 



3E(Y) 2(E(Y)) 2 J E(Y) J 2 



(note the added term in the coefficient of -^)- Therefore to exclude this special 
case, we also assume 

(A2) There exists some even k for which pk > 0. 

Now we add two technical assumptions. First, our proof of (|1.2p makes use of 
the assumption that 

(A3) There exists an N > such that pk = for k > N, 

that is, the jump when the origin is reached can only be to a finite number of 
sites. We note that this ensures that the moments E(Y), E(Y 2 ), E(Y 3 ) which 
appear in (|1.2I) will all be finite, but it is natural to conjecture that (|1.2p will 
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remain valid even when there are an infinite number of positive Pfc's, as long 
as these moments are finite. One may hope to prove (jl.2l) in the more general 
case by approximating an infinite sequence of pk's by finite ones, but some 
estimates will be needed in order to control the error term O ( ) during this 
approximation process, and we do not address this issue here. 

To formulate our second technical assumption, we introduce the generating func- 
tion for the probabilities pk ■ 

oo 

h(x) = j2p kxk > ( L4 ) 

fc=0 

and the associated function 

4>{x) = x 2 + 1 - 2xh(x), (1.5) 

which will play an important role. By (jl.ip . the series in fj 1 .4[) converges at least 
for \x\ < 1, so that h, <p are holomorphic in the unit disk. Under the assumption 
(A3), h and <f> are in fact polynomials. We will use the assumption 

(A4) The polynomial <f> does not have repeated roots. 

Note that this is a 'genericity' assumption, in the sense that for given degree, 
a generic polynomial of that degree will satisfy it (the polynomials failing to 
satisfy it are those with discriminant 0). This assumption can be eliminated 
by approximating polynomials with repeated roots by generic polynomials with 
distinct roots, but this requires careful control of error estimates, which we shall 
not carry out here. However, we do believe that the removal of assumption (A4) 
is an easier matter than the removal of assumption (A3). 

We can now state our main theorems 

Theorem 1 Under assumptions (Ai)-(A^), we have as n — > 00. 

Theorem 2 Under assumptions (A3), (A±) and 
(Ag) Pk = for all even k, 
we have il.S\) as n — > 00 . 

The proof of these theorems proceeds by the following steps 

• Generating functions encoding the expected values are derived (section^. 
In this section we also explore the connection to the spectral analysis of 
the transition operator associated to the process. 

• In order to study the asymptotics of the coefficients in the power series 
expansion of the generating functions, these functions are decomposed into 
sums of 'simple' components (section [3]). 

• Asymptotics of the coefficients of the power series of the 'simple' compo- 
nents are obtained using classical results on the relations between singu- 
larities of analytic functions and the behavior of their Taylor coefficients 
(section rj| . 
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• Everything is combined to obtain Theorems [T] and [5] (section [5]). 

We note that the generating functions for the probabilities, which easily imply 
the generating functions for the expected values, have already been derived in 
[5]. However, we give a complete derivation of these generating function in sec- 
tion [2j because we use a different approach, which relies on functional-analytic 
rather than probabilistic concepts. This approach reveals that the zeros of the 
function (f> defined by f) 1 .5[) are related to the eigenvalues and resonances of the 
operator transition operator L associated to our problem, and this fact puts the 
use of these zeros in the decomposition of the generating function into simple 
components, which a key element of the derivation of the asymptotic formulas, 
into a more general perspective. Thus, while our final results, Theorems 1 1121 
contain no mention of eigenvalues or resonances, and the use of the terminol- 
ogy of spectral theory is not essential for the proof of the of these result, the 
connection with spectral concepts enriches our understanding. 

Several directions for further work are suggested by the results and techniques 
presented here: 

• It will be of interest to eliminate the technical assumptions (A3),(A4) 
imposed here. 

• The asymptotic formulas are derived here by a technical calculation. It 
would be interesting to obtain a more intuitive understanding that explains 
the particular form that the terms in the asymptotic formula take. For 
example, is there a way interpret the constant shift term E ^(l^ 1 ? 

• The decomposition technique developed here should be of use in treating 
other stochastic processes in which a simple process (such as the simple 
random walk) is locally perturbed (in our example by the jump at the 
origin) . 

• While the relations between random walks (and other Markov processes) 
and eigenvalues is a huge topic (e.g. EJ [7]), it seems that the role of 
resonances in relation to stochastic processes has not received much at- 
tention. Resonances are most familiar to those studying quantum physics 
[5]. It would be interesting to further explore the role of resonances in 
stochastic processes. 

2 Derivation of the generating functions 

In this section we derive the explicit expressions for the generating functions of 
the expected values which are of interest to us. The results in this section are 
valid for any sequence of probabilities pk , and we do not need the assumptions 
(Ai)-(A4) introduced above. 

For each j > 0, we define the generating function 




(2.1) 



n=0 
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We will prove that 



Lemma 1 

HAz) = 2E(Y) P \ ' xx + j , 

JV ' V ' 1-z 4>{p{z)) J l-z' 

where 



p{z) = V - . (2.2) 



id <j) is defined by il.5\) . 



Let us note that throughout this paper, the function \fx refers to the branch of 
the square root defined on C \ (— oo, 0) with \/T = f . 

Lcmma[T]will follow from an expression for the two- variable generating functions 
encoding the probabilities 



oo oo 



F J (w,z) = J2Y, P ( X " = k \ X °=ti wkzn > 3>0- (2-3) 

Lemma 2 

F(u: z] _ 2 (p(z)P+VM~ ( Hp(z)V+ 1 
j[,Z) ~cf>(p(z))' z(w 2 + l)-2w ' [ > 

We note that equivalent expressions are already proved in [S]. 

To derive Lemma [T] from Lemma [2j we note that, by differentiating (12.31) with 
respect to w, and then setting w = 1, we obtain 

„ oo oo 

^J-(M) = 52z n *£p(x n = k\x = j)k 

j=0 k=0 
oo 

= Y / E(X n \X =j)z n = H j (z). (2.5) 

n=0 

Differentiating (|2.4p with respect to w, we have 

dF^ 2 (p(z)y+ 1 ^H-(j + l)^(z)V 

dw (W,Z> ch(p(z)) z(w 2 + l)-2w 

2 [(p(z)) j+1 <j>(w) - (f>{p(z))wi +1 ][2wz - 2] 
~ <j>{p(z)) [z{w 2 + 1) - 2w] 2 

and substituting w = 1, taking into account that 

= 2x - 2x/i(:c) - 2h'(x), 

and that 

oo oo 

h(l)=J2Pk = h h'(l) = J2kPk = E(Y), 

k=0 k=0 
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so that 

4>(1) = 2 - 2h(l) = 0, <f/(l) = -2h'(l) = -2E(Y), (2.6) 
we get, using (|2.5p 



9w ' 1 — 2 (f>(p(z)) 1 — z' 

which is the result of Lemma [TJ 

Let us note here that Lemma [T] immediately gives the 

Proof of Proposition]]} If po — 1 then = 0, so Lemma Q] gives 

1 



that is, in view of the definition (|2.ip . E(X n \ Xq = j) = j, for all j, n. ■ 

We now give the 

Proof of Lemma We introduce the functions 

oo 
fe=0 

so that Fj(w, z) defined by (|2.3p can be written as 

oo 

F j (w,z) = J2hn(™> n - (2-7) 

It will be useful to consider the functions fj. n (w) as elements of the Hilbert space 
H 2 of functions analytic in the unit disk with square-summable coefficients, with 
the inner product 

We have 

f ji0 {w)=w j , (2.8) 
and, in view of (|1.2j) . we can express fj, n +i i n terms of fj n through 

fi, n+l = L[f jtn ] (2.9) 

where L : H 2 H 2 is the bounded linear operator defined through its action 
on the basis {iy J }^. of H 2 by 

L[w j ] = ^w j - 1 + ^+\ j>l, 

OO 

L[l] = h(w)=Y J PkW k . 

k=0 

Let us note here that fj n encodes the state of the system started at j at time 
n, since the coefficients of f^ n give probability of being in each place at that 
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time. By (|2.9p . the operator L describes the dynamics of the system, that is 
the transition from the state at time n to the state at time n + 1. L may 
then be called the transition operator, and it is quite natural that, as we shall 
see below, the spectral analysis of L is closely related to the behavior of the 
stochastic process. 

The utility of formulating the operator in the Hilbert space H 2 is that L can 
be represented in a simple form as 

L[f] = + w)[f(w) - /(0)] + f(0)h(w). (2.10) 

By ([2J]) and fl2jj|) we have 

f j<n (w) = L n {f 3 , ] = L n [w j ]. (2.11) 
If \z\ < \\LW~ 1 then the resolvent 

R(z) = [I-zL]- 1 
is well defined, and we have the Neumann expansion 

oo 

71=0 



In view of (|2.11[) we have 

oo 

Fj{w, z) = J2 z n L n [m>] = R{z)[w j }. (2.12) 

n=0 

Therefore the function Fj(w,z) can be computed by computing the resolvent 
[I — zL]~ l , which we will now do. 

Let us then assume that /, g e H 2 with g ^ and 

[I - zL]f = g 

that is, by ([2~TU|) . 

f(w) - z^w- 1 + w)[f(w) - /(0)] - zf(0)h(w) = g(w). 

rearranging, we have 

2wg(w) + zf(0)[2wh(w) - (w 2 + 1)] zf(0)<j>(w) - 2wg(w) 



/H 



z(w 2 + 1) - 2w z(w 2 + 1) - 2w ' 

(2.13) 

where <f> is given by (|1.5I) . At this stage we have still not determined f(w) 
completely, since the right-hand side of (I2.13[) contains /(0), which still needs 
to be determined. To do so, we observe that in order for the function / defined 
by (|2.13[) to be an element of if 2 , it is necessary that it will not have singularities 
in the unit disk. We note that the quadratic equation 

z(w 2 + 1) - 2w = 
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has two roots, given by 



w 1 =p(z) = , w 2 



l + Vl^ 



: p(z) 
and for \z\ < 1 it is easy to check that 

\p(z)\ < 1. 

We may therefore rewrite (|2.13p in the form 

zf(0)(p{w) - 2wg{w) 



fH 



Z[w-p(z)][w- ^y] 



(2.14) 



(2.15) 



(2.16) 



Since, by (|2.15l) . w = is outside the unit disk, the only potential singularity 
is at w = p(z), and in order for this not to be a singularity we must have that 
the numerator vanishes at w = p{z), that is we must require 



zf(0)<f>(p(z)) - 2p(z)g(p(z)) = 0. 
Solved for /(0), (I2T71) gives 

2p(z)g(p(z)) 

Substituting (|2.18p back into (|2.13l) we obtain 

2 P{z)g(p{z))(j){w) - 4>(p(z))wg(w) 



R{z) 



fW 



z(w 2 + l)-2w 



(2.17) 



(2.18) 



(2.19) 



The function f(w) defined the right hand side of (|2.19p is defined and holomor- 
phic in P whenever 

<j>(p(z)) ? 0, (2.20) 

and, since p(0) = and <fi(0) = 1, this condition holds for \z\ sufficiently small. 
If, in addition to (I2.20p . we assume that 



z £ C\ ((-oo,-l] U [l,oo)), 



(2.21) 



then we can check that / £ H 2 , by using the integral form of the H 2 norm: 



\h 2 



1 

2Wo 



1 



1 



< 



47T \^{p{z))\ 2 
1 1 



\f(e l8 )\ 2 de (2.22) 
>T p(z)g(p(z))<P(e ie ) - <t>(p(z))e ie g(e ie ) 2 



cos(9)z — 1 



dO 



1 



47T \4>(p(z))\ 2 ee[o,2^] | cos(9)z - 1| 



p(z)g(p(z))^ e ) - J>(p(z)y e g(e ie ) 



d9, 



the maximum term being finite because of (12.211) and the integral being finite 
because of the assumption that g £ H 2 , and the fact that <j) £ H 2 , which follows 
from (|l.ip . This estimate also shows that R{z) is bounded. 
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In particular, taking g(w) = w j in (|2.19l) . and using (I2.12[) , we get 



Fj{w,z) = R{z)[w J \ = . 2.23 

(p(p(z)) z{w z + 1) — 2w 



We now make some comments regarding the spectral analysis of the transition 
operator L introduced above. We have seen above that the resolvent R(z) = 
[I — zL]^ 1 exists when z satisfies (|2.20p . (|2.21l) . Since the spectrum a(L) of L is 
the set of values of - such that R(z) does not exist, we have thus shown that 



<7(L)c{i|#p(z)) = 0}u[-l,l]. 



. z 

In fact we can see that the above containment is an equality: 
(i) The values - for which 



z e 



:\((-oo,-l]U[l,oo)), <f>(p(z)) = 0, (2.24) 



are discrete eigenvalues of L: indeed, assuming (|2.24[) and setting g = in 
(|2.16[) . we have that, for any choice of /(0), the function 



(W-P(Z))(W- ^jy) 

is holomorphic in B (note that <f>(p(z)) — implies that the apparent singularity 
at w = p(z) is removable), / € H 2 by a calculation similar to that in (I2.22p 
(using the assumption z £ (-co, —1] U [1, oo)), and [I — zL]f = 0. 

(ii) The shortest way to see that [—1,1] C <r(L) is to note that L is a compact 
(indeed a rank-one) perturbation of the discrete Laplacian, which is well known 
to have the essential spectrum [—1,1], so that by Weyl's theorem the essential 
spectrum of L is also [—1,1]. More directly, when z £ (— oo, — 1] U [l,oo) the 
function / defined by (|2.19l) will have a pole of order 1 on the boundary <9B, 
which will prevent / from belonging to H 2 . 

We thus have the following decomposition of the spectrum of L into the essential 
and discrete spectra: 

(x(L) = (r e3S (L)Ua diac (L) (2.25) 
a ess (L) = [-1,1] 

<JMsc(L) = {- I « eC\ (-oo.-l] U [l,oo), p(z) G <f>~\0)} 
The function p is a conformal mapping 

p:C\ ((-oo,-l]U[l,oo)) ->B, (2.26) 

its inverse given by 

p- 1 (o) = ^—, «£». (2.27) 
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We thus see that those zeros of <j> which are in the unit disc are in one-to-one 
correspondence with the eigenvalues of L. Wc furthermore claim that the zeros 
of <j) which are outside the unit disk correspond to the resonances of L (see [5] 
for an overview of resonances). To understand this, consider the behavior of 
the expression on the right-hand side of (|2.19l) . giving the resolvent of L, as 
z approaches the set (— oo, — 1) U (l,oo) from above or from below, that is as 
- approaches the essential spectrum of L. Since the expression 1 — z 2 then 
approaches the segment (— oo,0), where the square root has a branch cut, the 
limits as z approaches (— oo,— 1) U (l,oo) from above or from below will be 
different, but in each case analytic continuation is possible. Indeed, since the 
analytic continuation of y/x as x crosses the negative axis is —\/x, the analytic 
continuation of p(z) as z crosses the set (— oo,— 1) U (1, oo) will be (see 
(|2.14[> ). so that the expression 



1 - KikdMv) 



R(z)[g] - f(w) - —j- z(w 2 + 1) _ 2w ( 2 - 28 ) 



represents the meromorphic continuation of R(z) across (— oo,— 1) U (l,oo) to 
the 'second sheet'. Of course R(z) is no longer a bounded operator from H 2 
to itself, indeed the function f(w) has a singularity at w = p(z) 6 D. The 
important point for us is that R(z) has a pole when <^(^y) = 0, and since 

maps C\ ({— oo, —1] U [1, oo)^ to the exterior of the unit disk, we have that the 



set 



res{L) = {~ | J- € ^(0)} 



is now identified as the set of resonances of L, according to the definition of 
resonances as the poles of the meromorphic continuation of the resolvent. Wc 
note that when p(z) effln (/>~ 1 (0), then \ e [-1,1] is a resonance embedded 
in the essential spectrum. In fact - = 1 is always such a resonance, since 
p(l) = 1 G (/> _1 (0) - see Lemma [5] below. 

The zeros of <p will play an important role in the analysis carried out in the 
following sections, aimed at obtaining the asymptotic formulas given by The- 
orems 11121 although these zeros do not appear in the final formulas. The fact 
that these zeros corresponds to the eigenvalues and resonances of the operator 
L, as explained above, thus gives an interesting perspective. 

We shall now show that: 

<Jdi S c(L) C D. 
In view of (I2.25[) . this is equivalent to statement 

Lemma 3 

0(p(*)) = o => M>i. 

Proof: Assume that 4>{p{z*)) = 0. Set a* = p{z*). By (jH^SD we have \a*\ < 1, 
and by (|2T2"71) we have 

2a* 



a* 2 



z*. (2.29) 



11 



Since <j>(ce*) = we have, by (|1.5[) . that 



a* 2 + 1 

ft («*) = ^~- (2-30) 



Combining (|2~2l?|) and (g3D| we have /i(a*) = p-, hence 
Since \a*\ < 1, we have 



(2.31) 



|h(a*)| = I 5>*a**| ^ E^l a *! fc < X> = L 



' k=0 k=0 k=0 

Together with (12.311) . we thus have \z*\ > 1, as we wished to prove. ■ 
An important immediate consequence of Lemma [3] is 

Lemma 4 The functions Hj(z) are holomorphic in D. 

The following simple result will be used in later developments: 

Lemma 5 (i) x = 1 is a zero of <p(x), and if (A\) holds then it has multiplicity 
1. 

(ii) x — —1 is a zero of <j){x) if and only if (A\) holds. 
Proof : 

(i) Since h(l) = 1, we have 0(1) = 2 - 2h(l) = 0. Since <f/(x) = 2x - 2h(x) - 
2xh'(x), we have 

oo 

cf>'(l) = -2h'(l) = -2^ k Pk = -2E(Y). 

k=0 

If (Ax) holds then E(Y) > and thus tf/(l) < 0, so that the root x = 1 has 
multiplicity 1. 

(ii) If (Arj) holds then h(x) contains only odd powers of x, so h(— 1) = — 1, hence 
0(-l) = 0. Conversely, if 0(-l) = then ft(-l) = -1, that is 



£ Pfc (-i) fc = -i. 



fc=0 

By (jl.ip . this can hold only if p k =0 when fc is even, that is only when (A?,) 
holds. ■ 

Assuming, as we shall henceforth do, that (A3) holds, so that h(x), hence also 
4>(x), are polynomials, we have, denoting the degree of h (that is the maximal 
value of k for which pt- ^ 0) by N: 

• If N > 2 then the degree of (f> is N. The sum of the number of eigenvalues 
and resonances, counting multiplicities, is thus N. 
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• If N = 1 then <b{x) = x 2 + 1 — 2x[po + (1 — po)x], so the degree of 4> is 2, 
except in the special case 

V^=V\ = \ M^) = - + -x & <f>{x) =l-x. (2.32) 

In this special case the only resonance is z = 1, while otherwise there is 
an additional resonance given by z = po + 2p Q P -i • 

• If N = then = 1, that is po = 1, so that </>(x) = (x — l) 2 , so there is 
a double resonance at z = 1. This case is the one excluded by (^4i), and 
has already treated in proposition [TJ 



3 Decomposition of the generating functions 



Given a function f(z) holomorphic in a neighborhood of 0, we use the standard 
notation \z n ]f(z) to refer to the coefficient of z n in the power series correspond- 
ingto/, that is [*"]/(*) = £/ (n) (0). 

From Lemma Q] we have 

i p{zy +1 - 



E(X n | X = j) = [z n ][H 3 (z)] = 2E{Y) ■ [z r ' 



l-z cj){p{z)) 



(3.1) 



The problem we are faced with, then, is to estimate the power series coefficients 
on the right-hand side of (|3.1j) . 

By Lemma |Sfi) , the polynomial (f> can be decomposed as 

4>{x) = (x - 1)^(1), (3.2) 



where 4>(x) is a polynomial. Except in the special case given by (|2.32l) . the 
degree of ip is at least 1, so that it has at least one zero. Since we shall use the 
zeros of ip below, we will now exclude the special case (|2.32l) . and treat it by a 
separate calculation in section [51 which will show that Theorem [T] is still valid 
in this case. 

The first step is to perform a partial-fraction decomposition of the rational 
function ^^y- At this point we are using the assumption (Ai), so that the zeros 

of ip are non-degenerate, which allows us to decompose as 



tp(x) 



y -1 

^-^ ib' (a) x 



(3.3) 



and therefore 

x 



4>{x) 



E 



1 



1 



ib'{a) x — 1 x — a 



4>(a)=0 

so that, substituting x = p(z) 



E 

0(q)=O 



1 



1 



p(f) 

4>( P (z)) 



1 



ip'(a) a — 1 



a 



x — a 



1 1 



p(z) - 1 



(3.4) 
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Using the definition (|2.2[) of p(z), we have 
1 a 1 



az 



1 — z p(z) — a 1 — z 1 — vT" 



a 1 — az + \/\ — z 1 
1 - z [a 2 + l)z - 2a 



a 



[a -If 



(a 2 + l) 



{a 2 + l)z - 2a 1 - z \j 1- z 



1 — az + yl — z 2 1 — az 1 1 + z 



and 



l-z /£>(*) -1 1-z 1-vT 



z* — z 



1 Vl+I 



so 



1 



1 



p[z) — a p(z) — 1 



1- z 

a (a-1) 2 + (a 2 + a) z + VT^z 2 } la + 1 1 



(a -1)2 



so that (13.4[) gives 



(a 2 + l)z - 2a 



2a- 1 1 -z 



1 + z 1 yl + i a 



a-l) 2 V 1-2 2(1-2)1 



1 p(z) 



E ,w 



1 



1 - z 4>(p(z)) 

(a - l) 2 + (a 2 + 1)[(1 - a)z + y/l - z 2 \ 



i/)(a)=0 



i>'(a) L(a- l) s 



(a 2 + l)z - 2a 



a /l + 2 la + 1 1 



1 vm 



a 



(a - l) 3 V 1 - z 2 (a - l) 2 1 - z 2(a - 1) (i - z )i (a - l) 2 
We thus have the decomposition 



1 a(a 2 + 1) 



2 ^ E w7 



V>(a)=0 



V>'(a) (a - 1)' 



(3.5) 



where 



2E ^ E ^ 



i/>(a)=0 



j/>(a)=0 

^3 = W E ^ 

c 4 



V>'(a) (a-l) 3 ' 

1 a+ 1 

ip'Jaj {a-l) 2 ' 



1 1 



i/>(a)=0 



ij)'(a) a — 1 ' 



2i ^) E ^ 



i/>(a)=0 



<//(a) (a - 1) : 



(3.6) 
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and K a (z) is given by 



K ^ Z >~ (a 2 + l)z-2a • [6J) 

We now show that the numbers C\ , Ci , C3 can be explicitly evaluated as follows 
(C4 can be similarly evaluated, but it will not be needed): 



Lemma 6 Cx,C2,C 3 , defined by S3. 6\) are equal to 

1 E(Y 3 ) (E(Y 2 ) - l) 2 



C 2 

c 3 

Proof: Substituting x = 1 into (|3.3p we obtain 

1 1 



12 6E(Y) ME(Y)) 2 
E(Y 2 ) - 1 
2E(Y) ' 

1 

2' 



^ ib'(a)l-a ibilY ^ 3 ' 8) 



Cs = -E(Y)-±-. (3.9) 



>(1) 

Differentiating p. 31) we get 

J^iL = v — l — (310) 

(V>(*)) 2 ^J'(a)(x-a) 2 ' {6 - [J> 



and substituting x = 1, we get 



V -i ^ = ( 3ii) 

^ Q ^(a)(l-a)3 (V,(l)) 2 ' 1 • J 

From (13.81) . (|3.11[) and the decomposition rrz^~p = — T^U we nave 

C 2 = - J B(y) J] -i ° + \ =-2E(Y)-£^+E(Y)-L,. (3.12) 

Differentiating (|3.10[) we have 

fWW I fW = „ v 1 1 



and substituting x — 1 we get 

1 

^'(a) (q — l) 3 " 2(^(1)) 3 



V 1 1 _ ^(1)^(1)-2(^(1)) 2 



i/)(a)=0 
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From the decomposition ( a "^a = (a-i)' 1 + (a-i) A anc ^ P-HP - (QUE} we get 

Cl ^ 2i ^) E v^(T^F (3 - 14) 

f(l)^(l)-2(#)) 2 + 2^W(l) 

Differentiating (|3.2p we have 



${x) = yj(x) + (x-l)4>'(x), ())"{x)=2^{x) + {x-l)i)"(x), 
4>"'{x) = 3i>"(x) + (x-l)i)'"(x), (3.15) 

and setting x — 1 in (13. 15[) we have 

V.(1)=0'(1), ^'(1) = ^"(1), ^"(1) = ~<T(1), (3-16) 
so that (l3J)) . ((3~T2|) . (l3TT4t can be rewritten as 

Ci = E(Y) 



(^(l))3 

Differentiating (jl.5p we have 

0'(V) = 2ai - 2h(x) - 2xh'(x), <j>"{x) = 2 - Ah'{x) - 2xh"(x), 
<j/"(x) = -6h"(x) -2xh"'{x), (3.18) 

and substituting x = 1 into (|3 . 1 8|> . noting that 

N N 



fc=0 fe=0 



A' 



h"(\) = *(* - i)p* = £(^ 2 ) - ^(n 

fe=0 

N 

h"'(l) = jfe(Jfe - l)(fc - 2) Pk = E(Y 3 ) - 3E(Y 2 ) + 2E(Y), 



fc=0 

we have 



4>'{l) = -2E(Y), <P"(1) = 2(1 - E(Y) - E(Y 2 )), 

fi"(l) = 2{E{Y)-E(Y 3 )). (3.19) 



Substituting (|3.19p into (|3.17[) . we get the results. 
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4 Asymptotics of some power series coefficients 



In this section we derive some results on the asymptotics of the coefficients 
of power series expansions of the basic components which were obtained by 
decomposing the generating function Hq(z) in section [31 

A prototype result is (see [4], Theorem VI. 1) 
Lemma 7 For ueI\N, 

cj(oj + 1) 



[z n \{i- z y = 



r(- w ) 



2n 



O 



We need to deal with more complicated functions which are holomorphic in the 
unit disk but have singularities of the type (1 — z) u on its boundary. 

Let / be a function holomorphic in the unit disk D. A point £ with |£| = 1 is 
called a singularity if there is no neighborhood of £ such that / can be continued 
analytically to this neighborhood. A singularity £ is called algebraic with weight 
lj, where lj £ R \ N, if there exist holomorphic functions A(z), B(z), defined in 
a neighborhood INT of £, so that 



z G 



'HN^ f{z)=(l--XA{z) + B(z) 



(4.1) 



The following is Darboux's Theorem (the form given here is taken from pQ): 

Lemma 8 Assume f is holomorphic in D and has only a finite number of sin- 
gularities on its boundary, all of them algebraic, and let the minimal weight of 
the singularities be lj. Denote the singularities with weight lj by £i, • • • , £ m) and 
let Ai(z), Bi(z) (1 < i < m) be holomorphic in a neighborhood o/£,; so that 
holds with £ = £j. Then, asn->oo, 



[*"]/(*) = 



1 



An immediate corollary of Lemma [S] is 

Lemma 9 If j, k are integers with j > k, and k is odd, then, as n 



(1 + zp{l - z)* 



2i 



1 



r(-|) 



n 2 



1 
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Proof: £ = 1 is the singular point with minimal weight, co 



A(z) = (l + z)i, 



B(z) — 0. Hence the result follows from Lemma [51 ■ 

We now derive the estimates that we need for obtaining our main results. 

Lemma 10 As n — > oo, 



1 



1 



_2 3_ + n 



/ = 

/ > 1. 
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Proof: Set 



/(*) 



1 + Z ,4 2\I 14 . t+1 . , . 1 

(1 - z 2 )^ = (1 + z)~(l - z) 

1 — z 



If Z > 1 is odd then this function is a polynomial, so that its coefficients vanish 
for large n, so the 0[ -V ) estimate is trivially true. 

We now consider the case that I > 2 is even. Then Lemma [§J with j =1 + 1, 
k = I — 1 gives 



1 1 "t — n ~T~ \ n ~T~ / V «. 2 



1 



r(-— ) n 2 

There remains the case 1 = 0. We write 



1 

n 2 



1 

— 1 



where 



/(*) = 



0(2) 



1 + z 



1 - 2 
fc(z) 



= g(z) + V2^=L 



so that k(l) — 0, implying that 



yT - z 



(4.2) 



where 



</(*)= i4(*)(l-«) ! 

fc(z) 



A{z) = 



1 - z 



is holomorphic in a neighborhood of z = 1. Thus <? has an algebraic singularity 
with weight \ at z = 1. Clearly it also has an algebraic singularity with weight 
h at z = — 1 . From Lemma |8] we therefore have 



By Lemma [7] we have 



[z n ]g(z) = o( — 

v n2 



Vl - zt «2 V Vn 



(4.3) 



(4.4) 



Combining (|42 )) .([43 |) and (|4~4|) we obtain 



[z n ]f(z) = ^L -±= + 0(\ 
Jit Jn \ n 5 



Lemma 11 is n -> 00, 



[* n ] 



(1-z 2 )^ 



I = 



%* + °(j)' < = ' 

O(-V), ,>2. 
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Proof: Set 



/(*) = 



1 - z 



= (H-Z)3(1-Z)T 



When Z = we have 



When Z > 2 is even then /(z) is a polynomial of degree j, so the coefficients 
vanish for large n, and the 0[ -kr J estimate is trivially true. 

Assuming Z is odd, Lemma |H1 with j = Z, k = I — 2 gives 



[*"]/(*) 



2* 



r(i-i) 



In the case Z > 3, this gives the required estimate. 

If Z = 1 then f(z) = J jzr^, and the required estimate follows from the case 
Z = of LemmaHni ■ 



Lemma 12 As n — >• oo, 



{i- z y 



= < 



2 + oU-Y 

4 1 

V27T V" 

v V n 2 



Z = 
Z = 1 
Z = 2 
Z > 3. 



Proo/: Set 



/(*) = #^r(l -,¥ = (1 + ^(1-,)-. 

(1-Z)2 

If Z > 3 is odd then /(z) is a polynomial, so the coefficients vanish for large n, 
and the estimate holds trivially. 

If Z is even then from Lemma with j = 1 + 1, k = 1 — 3 it follows that 

1 



[*"]/(*) 



2^ 



r(- 



1-3 , _ 

— j n 2 



2 v n 2 ' 



Thus if Z > 4 then 



>]/(z) = 0(4r) 



If Z = 2 we decompose 



where 



fc(z) = (1 + z)* -2'. 
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and since fc(l) = we have 



g(z) = A{z)VT=z, A(z) 



1 - z 



where A(z) is holomorphic in a neighborhood of z = 1, so that z = 1 is an 



algebraic singularity of g with weight \ . 



with weight | . Therefore Lemma [5] implies that 

and using Lemma [7] we have 

[z n ]f(z) = [z n }g(z) + 2i[z n ] 
If 1 = 1 then 



T is also an algebraic singularity, 



1 



22 1 

/ 7T \Al 



o 



1 



so 



• M ; 1-2 

n/(z) = 2, »>i. 



If Z = we decompose 



where 



g( z ) = 



s/TT~z 

(l-z)i 

fc(*) 



V2 



\/2 1 



(I-z)^ 4 VT - T 



(4.5) 



, fc( z ) = VTTi-V2 + ^_(i- z ) 

(1 — Zj 2 4 



is chosen so that fc(l) = k'{\) = 0, hence setting 



g(z) = A(z)VT=z, A(z) 



Hz) 
(I-*) 2 ' 



we have that A(z) is holomorphic in a neighborhood of z = 1, so that 2 = 1 is an 
algebraic singularity of g with weight 5 . z = — 1 is also an algebraic singularity 
with weight i ■ Therefore lemma [5] implies that 



[z n ]g(z) = 0( — 
v n2 

Therefore, using (|4.5I) and Lemma [7] 

[z"]/(z) = [z n ]g(z) + V2[z n 
3 



1 



O 



1 

3" 

7? 2 



V2 



r( 



8n 



(l-z)t 
1 



V2 



n 2 



1 1 

2tt y/n 



O 



4 

V2 1 

Tf(|) 
1 

T 

n 2 



1 



V1^2 
l + O 



71/ / n 2 
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Lemma 13 For any I > 0, as n — > oo. 



v n2 / 



Proof: For Z even the function is a polynomial, so the result holds trivially. For 
I odd the result follows from Lemma [5] ■ 

Lemma 14 For a G C, Zei K a (z) be defined by \3. 7| ). 

fjj For anj/ a ^ ±1 wrai/i 0(a) = we have, for any integer I > 0, as n— » oo, 

For a = —1 we Ziawe, as n — > oo, 
[*»] [*•_,(*)(! -* 2 )4 



0(4)- 



o 



1 = 

i > i. 



Proof: Set: 



i (a-l) 2 + (a 2 + l)[(l-a)2 + VT^ 



/(*) = K a ( z) (l z^ = ^ ti ^ - ^ ^ — ^ _ z2 4 

(a z + l)z — 2a 

We first claim that K a (z), and hence also f{z), is holomorphic in D. Since 
the numerator and denominator are obviously holomorphic on this set, it is 
only necessary to show that the apparent singular point when the denominator 
vanishes, that is 

*o = — v (4-6) 
either satisfies \zo\ > 1, or is is removable. There are three cases to consider: 

(i) If \a\ < 1, then from (|2.27[) and (|4.6I) we have that p(zq) = a, and thus by 
the assumption that (j){ a ) = an d Lemma [3] we have that \zq\ > 1, so that the 
singularity is outside the unit disk. 

(ii) If \a\ = 1 then zq G (— oo, —1] U [1, oo), so the singularity is outside the unit 
disk. 

(iii) If \a\ > 1 then we will show that the numerator of K a (z) vanishes at z = zo, 
so that the apparent singularity at z = zq is removable. Indeed, |a| > 1 implies 
that 



so that 



I (a 2 - 1\2 



> o, 



1 



(note that the equality y x 2 — x holds only when $t(x) > 0). Therefore, substi- 
tuting z — zq into the numerator of K a (z) we obtain 



(a - l) 2 + (a 2 + 1) (1 - a)z + Jl - z 2 



21 



(a- l) 2 + (a 2 + 1) (1-a) 



2a 



i 2 - 1\2 



a 2 + 1 



= 0. 



We now assume that cn^il. 

The singularities of / on the boundary of the unit disk are £± = ±1, both 
algebraic singularities. 

If I is odd then the weight of both singularities is oj = | , since / can be repre- 
sented as 



/(z) = A + (z)(l - *)s + /(z) = A_(z)(l + z)s + B(z), 



where 



A ± (z) = 



[(a - l) 2 + (a 2 +!)(!- a)z](l ± *)i 



B(z) = 



(a 2 + 1)2 - 2a 



a +1 ., 2 \ '+ 



-(1-z 2 )- 



(a 2 + l)z - 2a 

are holomorphic in neighborhoods of ±1, respectively. Note that here we use 
the assumption a/±l, and that to conclude that B is holomorphic in neigh- 
borhoods of z = ±1 we use the assumption that I is odd. 

Lemma [8] thus gives 

[z n ]f(z) = C^— + o(^—) = 0(4r) ■ 

If Z is even then the singularities £± = ±1 have weights uj = ^jp, since / can be 
represented as 

f(z) = A + (z)(l - z) 1 ^ + B(z), f(z) = A_(z)(l + z)^ + B(z), 



A ± (z) 



(a 2 + l)(l±z) — 
(a 2 + l)z - 2a ' 



B(z) = (^-l) 2 + (^ + l)(l-^ (1 _ ^ 



(a 2 + l)z - 2a 



Therefore lemma [5] gives 



[z n ]f(z) = C-^ + of-^r) = o(4r 

There remains the case a = — 1. In this case we have 
/(z)=[2+ A /IHIl(l-z 2 )i 

Using Lemma [TOl we have 



1 + z 



1 - z 



4 /! + (-*) 



iri(i-^J = ^][a-(-) J ) l Vi^=4 
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= (-!)"[*"] 



1-z 
1 + z 



(1-z 2 )^ 



From Lemma [7] we have, for I > 1 odd, 



I > 1. 



1 



r(-|) 



1 + - 

n 



Of 



and the same estimate trivially holds for even I because the function is then a 
polynomial. 

Therefore if I > 1 we have \z n ]f(z) = O(-V), while for Z = we get 



[*»]/(*) = (-1)' 



1 



2tt Vn 







1 



5 Proof of the main theorems 



We now have the tools needed to finish the proofs of Theorems [T] and [21 

We consider first the case j — of Theorems 11121 From A3. 1 1) . (|3.5[) we have, for 
n > 1: 



E(X n | X = 0) 



[/][Io(z)] = d[/ 
2 ^ E ^ 



1 + z 



1 - z 



^'(a) (a - l) 3 



+ C 2 + C 3 [z n ] 
z n ][K a {z)\ 



(1-z) 
(5.1 



ip(a)=0 

(for n = one needs to add the term C4). 

For the first three terms in (I5.ip we have, from (|3.5[) and lemmas [TU1 and [T2l 



1 + z 



l-z 



= Cy 



-^ + c 2 + c 3 ( 



C 2 + C 3 [z r 
4 



Vl^K 



(l-z) 
1 1 



(5.2) 



2-7T 

2 _ E(Y 2 ) - 1 

fn A : '- 

2E(Y) 



/2?r \[n 



O 



2k 
1 



1 £(+ 3 ) , (E(Y 2 )-l) 2 ]J_ +0 (J_ 



3E(Y) 2(E(Y)) 2 J y/n 



3 



The behavior of the last term in (15.11) is determined by Lemma |T4j If (A 2 ) 
holds, then by Lemma [5jii) we have 1) 7^ 0, so the K-\(z) is not among 
the summands, hence, by Lemma [TW i) we have 



2E ^ E 777 w ^ [* B ][*«(*)] = o(4r)- 
_ V (a) (a — 1) Vri!'' 



i/>(a)=0 



(5.3) 
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If (A%) holds, then K—\(z) is among the summands, and by Lemma [T4Tii) it 
makes a higher order contribution so that 



ip{a)= 

(5.4) 

Note that = implies ft(-l) = -1, and (from (A$)) 



k—0 k odd 



SO 



[AD ^(-i) = -i^(-i) = -i[_2-2ft(-l)+2/ l '(-l)] = £(F), (5.5) 
hence (|5.4p gives 

2E(Y) Y - --. — [ z n ][K a {z)] = ^ + o(-A (5.6) 

ip{a)— 



From ([5TT |l .([53 1l and ([Q]l . wc get the case j = of Theorem [TJ From ([STTjl . flOJ) 
and (|5.6p , we get the case j = of Theorem [5J 

We now move on to the case j > 1. We note that, by Lemma [TJ we can write 

= p(^£r (*)+i T ^- > 



so that 

E(X n | X = j) - [z"]^(z) = [z n Mz)tH (z)]+j. (5.7) 

We shall find it easier to first estimate the coefficients [z n ] [z^ p(zY Ho(z)], and 
then use the obvious relation 

[z n MzyH (z)] = [z n +^p(zYH (z)}. (5.8) 



By the Binomial Theorem, 



zi P ( z y = [i - VT^y = £ ( J ) (-i)'(i - z 2 )K 

1=0 ^ ' 
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hence, using (|3.5p . 

[/][^w( 2 )]=c,t(;|(-i)'n 



2=0 



l + z 
1-z 



(l-z 2 )^ (5.9) 



2=0 

3 



(1-z 2 )^ 



1-z 



2=0 

3 



(-i)'M 



(1-z)' 



-(1-z 2 )^ 



2=0 



2£ ^) E 



^ P (-i)'[/]k( z )(i- 



V>'(a) (a- l) 3 W 

Examining the first sum in (|5.9p , we see from Lemma [TU] that the only term in 
this sum making a contribution up to order O ( -V ) is the term corresponding 
to / = 0, giving 



3 
J 

2=0 

2 1 



1 + 2 
1-Z 



(1-z 2 ). 



l + z 



1 - z 



2?r V" 



Of — 
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(5.10) 



Similarly, for the second sum in (I5.9P we see from Lemma[TT]that the only terms 
making a contribution up to order O ( -V ] are those corresponding to Z = and 



/ = 1, giving 



2=0 



£(;|j(-i)V][(i-* 2 )* r 

.(1-z 2 )^ 



(5.11) 



1 - z 



O-r =i-j 



2 1 



2tt \Jn V n 5/ 



For the third sum in (15.91) , we see from Lemma [12] that the terms making a 
contribution are those corresponding to I = 0, 1, 2, giving 



2=0 



i vT+z 



(l-z)t 



(l-z)t 



■ 3[z 



-(l-z 2 ). 



1 •/ ■ ,u m r VTT 



2tt 



(1-z) 

j' + (|+iCi-i)) 



(1-z) 

F (i-^ 2 )l+o(4) 

2 J V n, 2 i 



1 1 

2tt 



1 

71 2 



(5.12) 



25 



For the fourth sum in (|5.9p , Lemma Q2] implies that none of the terms make a 
contribution up to order O ( -\ ] , that is 



2=0 



o 



(5.13) 



The behavior of the last sum in (|5.9p is determined by Lemma [TU In case (A2) 
holds, so that a = — 1 is not a zero of ip (Lemma EJii)) we have from Lemma 

[Tjn) that all terms in this sum are of order O ( -\ ) , so that 



i/>(a)=0 



1 
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In the case that {A? 2 ) holds, we have the term corresponding to a = —1 in 
the sum, and in this case Lemma [TIT ii) tells us that the term corresponding to 
a = — 1 and I = gives us a contribution, so that 



^,(a)=o Y y ' v ; 2=0 v 7 



— 7. 2 



1 1 



8^'(-l) 



[*»][#_!(*)] + 



^ 1 (-1) 



n+l 



o 



(5.15) 

\r>. 2 / 



where in the last step we have also used (|5.5I) . 

Summing all these estimates, we have that, in the case (A 2 ) holds, (|5.9p - (|5.14p 
imply 



[z n ][z>p(zYH (z)] = ~ 7 =V^+ (C 2 - j) 

V 27T 



(5.16) 



+ ( Ci - jC 2 + - 



1 , iO'-ih 2 1 



2-7T V™ 



O 



3 / ) 

n2 ■ 



while in the case that (Afj) holds, (|5T9 ]) - ([5TT2) l and ([535)1 imply 



[* n ][^p(*yflo(*)] - (Ca - j) 

V Z7T 



(5.17) 



Assuming now that (A 2 ) holds, from (|5.8I) and (I5.16P we have 



2 1 

27r \fn 



01 



V Z7T 



(5.18) 



+ (C!-iC 2 + 2 
and in view of 



1 , i(i - 1; 



1 



o 



1 



/27T y/n + j \{n + j) 



3 1 



Vn + i = Vn\/l + - = Vn+^-7=+0(— ), 



1 



1 1 



1 + f 



2 

-i + o(4r) 



(5.19) 



2G 



([5TT5j) implies 



[z n ][p(zyH (z)] = JLy/K+tQt-j) (5.20) 
v 2ir 



So that together with (|5.7p we have 



| X = j) = -L + C 2 (5.21) 

V27T 



/ 1 i 2 \ 2 1 /I 

V 4 2/ V27T V n v n5 

and substituting the values of C\ , C% , C3 as calculated in Lemma [51 we obtain 
the result of Theorem [T] 

In the case that (A2) holds, we carry out the same procedure, starting with 
(|5.17[) . to obtain the result of Theorem [2J 



6 The case Po = Pi = -k 



It remains to treat the special case given by (|2.32l) . Since in this case <fi(x) = 
1 — x, Lemma [T] gives 



1 — z 1 — p(z) 1 — z 



VTTz" 

(l-z)t 



1 



l-z 



1 - v 7 ! -z< 



j'+i 



+ 3 



z j+l 
1 - Z 



2 



(1-2)4 



3+1 

E 

Z=0 



.7 + 1 



(-l)*(l-* 2 )*+j 



1 - x 



So that, calculating as in section [SJ we have that for n > j + 1 

3+1 



[z n ][z> +1 Hj{z)] = lY, 



3 + 1 
I 



(-!)'[* 



'r.n,r 7 ^±^ (1 _ za) * 
(1 - z) 



1 3+1 

Z=0 



i + i 



(-I)V' 



(1-z 2 )^ 



1 - 2 



11 / . . , 2 1 

-2(j + i)+j(j + iy 



2ir y/n 



2?r \fn. 



1 - (3 + 1) 
- 1 



2tt 



+ r-j 



2 M 

— ~r= +3 + 
2vr 

1\ 1 1 



2tt Vn 



1 

1 

7) 2 

fO 
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Therefore, using (|5.19l) . 



E(X n | X = j) = [z n ][H 3 (z)] = [z n+ ^][z^H 3 {z)] (6.1) 

2 , 1 / , h 1 1 

—=Jn +7 + 1 h 2r - - — — , 

+ b -t;)-s=-^ + 0| 



2/^-1, 
7^(^ + 2 ( - ? + 1) - 



n 2 



2 ^ 1 /. 2 . 1\ 1 1 



Since in the case under consideration we have E(Y) = E{Y 2 ) = E(Y 3 ) = ~, it 
can be checked that the result (I6.1[) is exactly that given by Theorem [TJ 
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